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I Abstract. Hadron spectra in Au+Au colhsions at RHIC are calculated by 

' hydrodynamics with initial conditions from the Color Glass Condensate (CGC). 

Minijet components with parton energy loss in medium are also taken into account by 
using parton density obtained from hydrodynamical simulations. We found that CGC 
provides a good initial condition for hydrodynamics in Au+Au collisions at RHIC. 
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^ ' High energy heavy ion colhsions involve different aspects according to the relevant 

^ ■ energy or time scale. There already exist many theoretical approaches to understand 
O ' numerous RHIC data. In this work, we consider in particular physics of gluon saturation, 
2 ' hydrodynamic evolution, and the energy loss of hard partons in the medium. Our goals 
^ ■ are to combine them and to take a step to a unified understanding of the dynamical 
aspect of high energy heavy ion collisions. 

At soft region < 2 GeV/c, where bulk dynamics is governed, hydrodynamical 
■ description P is successful in describing the elliptic flow at low pr, up to semi-central 
^ ■ collisions, and mid-rapidity j2] at RHIC. This is one of the strongest indications of 
an early thermalization of the quark gluon plasma (QGP) at RHIC. Hydrodynamics 
also predicts that the scaled elliptic flow, which is defined as the second harmonics 
V2 divided by initial spatial eccentricity e, becomes almost constant around 0.2 [3]. 
The experimental data reaches the hydrodynamic limit for the first time in central and 
semi-central collisions at RHIC energies |3]. On the other hand, minijets go across the 
expanding matter and lose their energies (jet quenching) in heavy ion collisions [S]. 
Observed large suppression of hadron spectra jE] and disappearance of the away-side 
peak in azimuthal correlation functions at mid-rapidity [7] have been interpreted as 
a consequence of parton energy loss in the medium. Indeed, Cronin enhancement of 
the hadron spectra and existence of the back-to-back correlation at mid-rapidity in dA 
collisions at RHIC [Hj support the importance of the strong final state interaction in 
AA collisions. The current RHIC data strongly suggest that the initial parton density 
is large. What is an origin of the large density in Au-|-Au collisions at RHIC? The bulk 
particle production is dominated by the small x modes in the nuclear wave function, 
where x is a momentum fraction of the incident particles. It is well known that gluon 
density increases rapidly with decreasing x by the BFKL cascade until gluons begin 
to overlap in phase space where nonlinear interaction becomes important These 
gluons form the Colour Glass Condensate (CGC) fUl- Remarkably, the CGC results 
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on the global observables such as the centrality, rapidity and the energy dependences 
of charged hadron multiplicities agree with the RHIC data ^1]. It has been shown 
that the classical wave function in the MV model contains Cronin enhancement and the 
quantum evolution in x makes the spectrum suppressed jT^j. 

From the above considerations, the CGC, hydrodynamics, and the energy loss of 
hard partons are key ingredients to describe the RHIC physics and must be closely 
related with each other. For example, the CGC could be a good initial condition for 
thermalization because it produces a large number of gluons. Thus these gluons are 
responsible to the large suppression of jet spectra. In this work we assume that the 
origin of thermalized partonic matter is the CGC in high energy heavy ion collisions, 
and use it as an initial condition in the hydro+jet model ^S]- With this approach we 
expect to get deeper understanding of the dynamical aspect of the heavy ion collisions. 
In addition, some of the problems which are inherent in a particular approach can be 
removed. We employ the /c^ factorized formula along the line of Kharzeev, Levin, and 
Nardi (KLN) ^T] for the computation of the gluon rapidity distribution which is given 



where Xi^2 = Pt exp(±y)/-ys with y and Pt are a rapidity and a transverse momentum 
of a produced gluon. We assume that the system of initially produced gluons reaches 
local thermailzed state at a short time scale. Although the produced gluons will reach 
the thermalized state through the dissipative processes in the realistic situations, the 
description of non-equilibrium phenomena is beyond the scope of the present paper. We 
also assume that the shape of the rapidity density distribution is not changed during 
the system is thermalized. Therefore, we take initial conditions from gluon distribution 
obtained from Eq. (^. Assuming Bjorken's ansatz y = f]s where ?7s is the space-time 
rapidity, we obtain the number or the energy density for gluons at each space-time point. 

The initial transverse energy per particle yields Ex/Ng ~ 1.6 GeV at ?/ = from 
Eq. (PJ. This is within a range estimated in a numerical simulation of the classical Yang- 
Mills equation ^3]. It should be noted that the assumption of the thermalization in 
CGC gluon distribution is to reduce the transverse energy per particle from Ex/Ng = 1.6 
to Ex/Ng ^ 1 within our present parameters. The effect of the hydrodynamic 
afterburner is to reduce the transverse energy per particles due to pdV work and 
yields {dEx / dy) / {dN / dy)\y=o = 0.54 GeV. The rapidity distribution becomes slightly 
wider. Our result supports that KLN calculation [TTj which is based on the assumption 
of parton-hadron duality is a good approximation on the rapidity distributions. In 
Fig. ^ pseudorapidity distributions of charged hadrons in Au + Au collisions at both 
y'sjvAf = 130 and 200 GeV are compared with the PHOBOS data ^3]. kr factorization 
approach in the CGC provides very good initial conditions for the hydrodynamical 
simulations which reproduce rapidity, centrality and energy dependences of multiplicity. 
It should be emphasized that it is very hard to find such a good initial condition which 
fits the data with the same quality as the CGC one presented here. 
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Figure 1. Pscudorapidity distributions of charged hadrons in Au + Au collisions at 
y/SNN = 130 and 200 GcV is compared to the PHOBOS data ^5]. Impact parameters 
which correspond to (A^part) from PHOBOS are used in the calculations. 




Figure 2. Transverse momentum spectra for negative pions (left) and nuclear 
modification factor Raa (right) in Au + Au collisions at y/sJjN = 200 are compared 
to the PHENIX data 



We now turn to the discussion of the high pt hadron spectrum. In our model, 
high Pt jets suffer interaction with the local parton density which is governed by 
hydrodynamic evolution. We only take into account parton energy loss in deconfined 
matter T > Tc. In Fig. |2l the result of the centrality dependence of the pion spectrum 
and the nuclear modification factor Raa integrated over pt > 4.5 GeV/c for |?7| < 0.35 
are compared to PHENIX data [12]. Our results only account for the data up to mid- 
central events and fail to reproduce data at peripheral collisions, because neither CGC 
nor hydrodynamics can be applied in the low density region. The centrality dependence 
is well described by assuming the number of participants scaling However, note 
that our density scales as p ~ — Trw:Q'i ~ — TrW\Ppa.Tt at mid-rapidity. By comparison, 
we also plot the LOpQCD calculations with the same initial condition assuming (1-|-1)D 
expansion of the system (p ~ 1/r). LOpQCD calculation without parton energy loss at 
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T < Tc = 170 MeV is consistent with our hydrodynamic result. However, if there is no 
restriction on the minimum temperature in the calculation of energy loss, we see that 
agreement with data becomes somewhat better. This may indicate the contribution 
from hadronic interactions at peripheral colhsions. 

In summary, CGC, hydrodynamics and the energy loss of hard jets have been 
integrated into one model and this dynamical approach describes the bulk properties of 
RHIC data. A systematic study within this dynamical approach is in progress. 
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